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Abstract 

Small diameter asymptotics is obtained for scattering solutions in a network of thin fibers. 
The asymptotics is expressed in terms of solutions of related problems on the limiting quantum 
graph F. We calculate the Lagrangian gluing conditions at vertices v 6 T for the problems on 
the limiting graph. If the frequency of the incident wave is above the bottom of the absolutely 
' continuous spectrum, the gluing conditions are formulated in terms of the scattering data for 

each individual junction of the network. 
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1 Formulation of the problem and statement of the results 

The paper concerns the asymptotic analysis of wave propagation through a system of wave guides 
when the thickness e of the wave guides is very small and the wave length is comparable to e. The 
problem is described by the stationary wave (Hclmholtz) equation 
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- e 2 Au = \u, i£f! E , (1) 

in a domain Sl e c R d , d > 2, with infinitely smooth boundary (for simplicity) which has the following 
structure: Q £ is a union of a finite number of cylinders Cj <e (which we shall call channels), 1 < j < N, 
of lengths lj with the diameters of cross-sections of order O (e) and domains Ji )£ , • • • , Jm,e (which 
we shall call junctions) connecting the channels into a network. It is assumed that the junctions have 
diameters of the same order 0(e). Let m channels have infinite length. We start the numeration 
of Cj iE with the infinite channels. So, lj = oo for 1 < j < m. The axes of the channels form edges 
Tj of the limiting (e — + 0) metric graph T. The vertices Vj £ V of the graph T correspond to the 
junctions Jj jS . 

The Helmholtz equation in f2 e must be complemented by the boundary conditions (BC) on <9fL. 
In some cases (for instance, when studying heat transport in f2 £ ) the Neumann BC is natural. In 
fact, the Neumann BC presents the simplest case due to the existence of a simple ground state (a 
constant) of the problem in fl £ . However, in many applications, the Dirichlet, Robin or impedance 
BC are more important. We shall consider (apart from a general discussion) only the Dirichlet BC, 
but all the arguments and results can be modified to be applied to the problem with other BC. 

An important class of domains fl e are self-similar domains with only one junction and all the 
channels being infinite. We will call them spider domains. Thus, if f2 e is a spider domain, then there 
exist a point x = x(e) and an e-independent domain f2 such that 
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Figure 1: An example of a domain f2 £ with four junctions, four unbounded channels and four 
bounded channels. 



e = {(x + ex) : x £ fi}. 



(2) 



Thus, fl £ is e-contraction of fl = 

For any D e , let Jj( v ),e be the junction which corresponds to a vertex v £ V of the limiting graph 
r. Consider a junction Jjf v ),s an d a U adjacent to Jji v ),s channels. If some of these channels have a 
finite length, we extend them to infinity. We assume that, for each v £ V, the resulting domain 0„ jE 
which consists of junction Jj( v ),e an d emanating from it semi-infinite channels is a spider domain 
(i.e., £l v , E is self-similar). This assumption can be weakened. For example, one can consider some 
type of "curved" channels, and the final results (with some changes) will remain valid. Simple 
equations on the limiting graph in this case will be replaced by more complicated equations with 
variable coefficients. However, even small deviation from the assumption on the self-similarity of 
f2„ e would make the statement of the results and the proofs much more technical. So, we consider 
only domains £l £ for which £l Vi£ , v £ V, are self-similar. 

Hence, the cross sections ojj >E of channels Cj. £ are e— homothety of bounded domains uij £ 
Let Xj.o < Xj t % < Aj.2--. be eigenvalues of the negative Laplacian — A^-i in ujj with the Dirichlet 
boundary condition on duij, and let {<fj, n } be the set of corresponding orthonormal eigenfunctions. 
The eigenvalues Xj jn coincide with the eigenvalues of — e 2 Ad-i in u>j, £ . In the presence of infinite 
channels, the spectrum of the operator — e 2 A in Q e with the Dirichlet boundary condition on dil £ 
has an absolutely continuous component which coincides with the semi-bounded interval [Ao,oo), 
where 



The equation ([T]) is considered under the assumption that A > A , when propagation of waves is 
possible. There are two very different cases: A — * A as e — > 0, i.e. the frequency is at the edge 
(or bottom) of the absolutely continuous spectrum, or A — * A > Ao, i.e. the frequency is above 
the bottom of the absolutely continuous spectrum. There are many results about the first case, the 
references will be given later. This paper concerns the asymptotic analysis of the scattering solutions 
for the Dirichlet problem in f2 e when A is close to A > Ao- 

If e — > 0, one can expect that the solution u £ of (p} in fl £ can be described in terms of the solution 
S = Ce(t) of a much simpler problem on the graph T. For example, if A^o < A < Aj,i for all j, then 
C satisfies the following equation on each edge of the graph 



min A 

l<j<m 



(3) 



e 2 d\ (t) 
dt 2 



(A-A,-oMt), 



(4) 
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where t is the length parameter on the edges. One has to add appropriate gluing conditions (GC) 
at the vertices v of T. These gluing conditions give basic information on the propagation of waves 
through the junctions. They define the solution ^ of the problem Q on the limiting graph. The 
ordinary differential equation (U|), the GC, and the solution depend on e. However, we shall often 
call the corresponding problem on the graph the limiting problem, since it enables one to find the 
main term of the asymptotics as e — > for the solution u = u e of the problem ([I]) in f2 E . 

One of the main difficulties in the problem under investigation was to find the GC, in particular, 
since the GC differ dramatically from those which were known in the case of A close to the bottom 
of the spectrum. 

Let us define the scattering solutions for the Dirichlet problem in fi e . We introduce local coordi- 
nates (t,y) in each channel C^ e with t axis parallel to the cylinder Cj jE , < t < L, and y £ i? n_1 
being Euclidean coordinates in the plane perpendicular to the t axis. The coordinate y is chosen in 
such a way that Uj iE = {(sy) : y E ujj E i?™ -1 }. For each j, the set {eT"^, n (^ )} is the orthonormal 
basis in L 2 (ojj i£ ) consisting of eigenfunctions of the operator — e 2 Ad-i- 

Let I be a bounded closed interval of the real axis which does not contain the points Xj, n , j < N. 
Thus, there exist rrij > f such that Xj, mj < A < \j, mj +i for all A E I. As will be seen from the 
definitions below, rrij + 1 is the number of waves which may propagate in each direction in the 
channel Cj_ e without loss of energy and with frequencies less than VA, A E I . We put rrij = — f, 
thus {Xj,m < n < rrij} is the empty set if A^o > A for A E I. 

Consider the non-homogeneous Dirichlet problem 

(-e 2 A-A)u = /, x E tt E ; u = ondQ E . (5) 

Definition 1 Let f E L 2 com {Q, e ) have a compact support, and A € I. A solution u of |5[) is called an 
outgoing solution if it has the following asymptotic behavior at infinity in each infinite channel Cj jE , 
1 < j < to: 

rrij 

u = Y / a 3 , n e l ^^ lt ^ n (y/e) + 0(e-^), 7 = 7 (e) > 0, (6) 

71=0 

Definition 2 A function ^> = 'I'^jj., f<s<TO, 0<fc< rrij, is called a solution of the scattering 
problem in O e if 

(-e 2 A - A)* = 0, x E fl £ ; * = on dfl £ , (7) 
and ^ has the following asymptotic behavior at infinity in each infinite channel Cj^ e , 1 < j < to: 

rrij 

= S s , j e- i >^ Kt <p s , k (y/e) + £ t j>n e i '^ Kt lf j, n ( y /s) + 0(e^% (8) 

n=0 

where 7 = 7(e) > 0, and <5 SJ is the Kronecker symbol, i.e. S s j = 1 if s = j, 5 s _j = if s ^ j. 

The first term in 1(5)) corresponds to the incident wave, and all other terms describe the trans- 
mitted waves. The incident wave depends on s and k, where s determines the channel, and s and k 
together determine the frequency of the incident wave. The transmission coefficients tj^ n also depend 
on s and k (i.e. on the choice of the incident wave), so sometimes we will denote them by tj' n . 

We introduce an order in the set of incident waves and corresponding scattering solutions and the 
same order in the set of transmitted waves. Namely, we number the incident waves in the channel 
Ci j£ taking them in the order of increase of absolute values of their frequencies, then we number all 
the solutions in the channel Ci^ e , and so on. With this order taken into account, the transmission 
coefficients for a particular scattering solution form a column vector with 



m 

M = J^irnj + 1) 



(9) 



entries. Together, they form an M x M scattering matrix 

T = {tfj, (10) 

where s, k define the column of T and j, n define the row. We denote by D the diagonal M x M 
matrix with elements -JX — Aj, n on the diagonal taken in the same order as above. The following 
statement can be useful in some applications, and will be proved in the next section (although it 
will not be used in this paper). 

Theorem 3 The matrix D 1 / 2 TD~ 1 / 2 is unitary and symmetric. 

The operator H — — e 2 A with the Dirichlet boundary conditions on <9f2 e is non-negative, and 
therefore the resolvent 

Rx = (s 2 A- X)- 1 : L 2 (Q S ) -> L 2 {Vl e ) (11) 

is analytic in the complex A plane outside the positive semi- axis A > 0. Hence, the operator R^t 
is analytic in k in the half plane Imfc > 0. We are going to consider an analytic extension of the 
operator R^i onto the real axis and in the lower half plane. Such an extension does not exist if R^i 
is considered as an operator in L 2 (Sl e ) since R& is an unbounded operator when A = k 2 belongs to 
the spectrum of the operator R\. However, one can extend R^ analytically if it is considered as an 
operator in the following spaces (with a smaller domain and a larger range): 

R k i : L 2 com (n e ) - Ll c (Q s ). (12) 

Theorem 4 (1) The spectrum of the operator H = — e 2 A in Cl e with the Dirichlet boundary con- 
ditions on dtt £ consists of the absolutely continuous component [Ao,oo) where Ao > is given by 
(0) and, possibly, a discrete set of positive eigenvalues {X J ' E } with the only possible limiting point at 
infinity. The multiplicity of the a.c. spectrum changes at points A = Xj, n , and at any point A, it is 
equal to the number of points Aj jTl , 1 < j < m, located below A. The eigenvalues A J,£ = X J for spider 
domains £l e do not depend on e. 

(2) The operator h!2\) admits a meromorphic extension from the upper half plane Imk > into 
lower half plane Imk < with the branch points at k = ± y/Xj, n of the second order and the real 
poles at k = ±V A s,e and, perhaps, at some of the branch points. The resolvent $12)) has a pole at 
k = ii/A^n if and only if the homogeneous problem with X — Xj_ n has a nontrivial solution u 
such that 

u= ^ a j,n<Pj,n(v/ e ) + ( e ~ 7 *)> xeCj.e, t^co, 1 < j < m. (13) 

j,n:Xj n =X 

(3) If f E L 2 com (VL £ ), andk = \/J is real and is not a pole or a branch point of the operator 
H12]) . and A > Ao, then the problem (0), is uniquely solvable and the outgoing solution u can be 
found as the Lf oc (fl £ ) limit 

u = R x+i of- (14) 

(4) There exist exactly M (see different scattering solutions for values of X > Ao such that 
k = VX is not a pole or a branch point of the operator U2\) . and the scattering solution is defined 
uniquely after the incident wave is chosen. 

Remarks. 1. Operator H = — e 2 A and its domain depend on e. One could use the term "family 
of operators" when referring to H. We prefer to drop the word "family", but one must always keep 
in mind that H depends on e. 

2. Existence of a pole of the operator (fT2")l at a branch point means that R^i has a pole at z = 
if this operator function is considered as a function of z = \/k 2 — Xj n . 
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3. One can not identify poles of the resolvent and eigenvalues of the operator based only on 
general theorems of functional analysis since we deal with the poles of the modified resolvent (fT!?|) 
which belong to the absolutely continuous spectrum of the operator. 

4. The eigenvalues A J:£ of the operator H can be embedded into the absolutely continuous 
spectrum, and can be located below the absolutely continuous spectrum. In particular, from the 
minimax principle it follows that H necessarily has a non-empty discrete spectrum below Ao if at 
least one of the junctions is wide enough. For example, non-empty discrete spectrum below Ao exists 
if a junction contains a ball B p of the radius p = re such that the negative Dirichlet Laplacian in 
the ball B r has an eigenvalue below Xq. 

Let us describe the asymptotic behavior of scattering solutions W = as e — > 0, A 6 I. Note 
that an arbitrary solution u of the equation ([1]) in a channel C JiE can be represented as a series with 
respect to the orthogonal basis {<Pj,n(y/£)} of the eigenfunctions of the Laplacian in the cross-section 
of Cj t£ . Thus it can be represented as a linear combination of the travelling waves 



e ±l » 3 Vj\n(y/e), 1 < n < rrij, 

and functions which grow or decay exponentially along the axis of Cj, e . The main term of small 
e asymptotics of scattering solutions contains only travelling waves, i.e. on each channel Cj. £ , any 
function \P has the form 

rrij 

* = *S = Ef^"^' + hne-'—^^nivh) + <k, (15) 

71=0 

where 

\r e s>k \ < Ce- 1 ^ 1 , 7>0, and d(t) = mm(i, lj - t). 

The constants a.j t n and f3j. n depend also on s, k and e. The formula (| 1 5|) can be written in a shorter 
form as follows 

rtij 
n=0 

where tpj — ipj(y/e) is the vector with components <pj tn (y/s), < n < rrij, and $j = <>j(t) is 
a {rrij + l)-vector whose components qj^ n are linear combinations of the corresponding oscillating 
exponents in t, i.e. <Tj satisfies the following equation: 

( £2 |i+^ = ' < £ < Z,, (16) 

where Dj is the diagonal matrix with elements yX— Aj.„, < n < rrij, on the diagonal. 

In order to complete the description of the main term of the asymptotic expansion (|15p . we 
need to provide the choice of constants in the representation of Cj,n as linear combinations of the 
exponents. Thus, 2(rrij + 1) constants must be chosen for each channel Cj^ £ . We consider the limiting 
graph r, whose edges Tj are the axes of the channels C Jl£ . Let <; be the vector valued function on T 
which is equal to qj on Tj . The vector ? has a different number of coordinates on different edges Tj 
of the graph T. We specify ? by imposing conditions at infinity and gluing conditions (GC) at each 
vertex v of the graph T. Let V — {v} be the set of vertices v of the limiting graph T. These vertices 
correspond to the junctions in £l £ . 

The conditions at infinity concern only the infinite channel Cj e , j < m. They depend on the 
choice of the incident wave and have the form: 

a. - { 1 if C?» = ( s ' fc ) l<i< m (17) 
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The GC at vertices v of the graph T are universal for all incident waves and depend on A. In 
order to state the GC at a vertex v , we choose the parametrization on F in such a way that t = 
at v for all edges adjacent to this particular vertex. The origin (t = 0) on all other edges can be 
chosen at any of the end points of the edge. Consider auxiliary scattering problems for the spider 
type domain fi„ iE formed by the individual junction, which corresponds to the vertex v, and all 
channels with an end at this junction, where the channels are extended to infinity if they have a 
finite length. We denote by T v the limiting graph which is defined by fl v . e . Definitions [IJ [2] and 
Theorem |4] remain valid for the domain Q v-£ . In particular, one can define the scattering matrix 
T = T v for the problem (TIJ) in the domain fi„, s . Let v%, V2,---vi, I — l(v), be indices of channels in 
O e which correspond to channels in f2„ i£ . Les us form a vector by writing the coordinates of all 
vectors <^ Us in one column, starting with coordinates of <; V1 , then coordinates of <^ V2 , and so on. Let 
us denote by D v (A) the diagonal matrix with the diagonal elements — X Vs ~ written in the same 
order as the coordinates of the vector q( v \ Let I v be the unit matrix of the same size as the size of 
the matrix D V (X). The GC at the vertex v has the form 

e[I v +T v ]D- 1 (X)j t &Ht) +l [I v -T v ]&\t) = 0, t = 0. (18) 

The GC (|18p has the following form in the coordinate representation. Let Z = Z(v) be the set of 
indices (j, n), where j are the indices of the edges of T ending at v and < n < rrij. Then 



E 

(i,n)£Z 



S j£ + tj'£{ v ) ( a ~ a j,«) 1/2 ^> + « 6j£-tj*(v) = at v, 

(s,k) € Z, 



where tj' n (v) are the transmission coefficients of the auxiliary problem in the spider domain £l Vy£ 
(i.e. are the elements of T v ), and 8*' n = 1 if (s, k) = (j, n), 5?' n = if (s, k) ^ (j, n). 

Definition 5 A family of subsets 1(e) of a bounded closed interval I C R 1 will be called thin if for 
any S > 0, there exist constants [3 > and ci, independent of S and e, and C2 = C2(S), such that 1(e) 
can be covered by c\ intervals of length S together with C2£ -1 intervals of length cie~^l e . Note that 
\l(e)\ -> as e -► 0. 

Theorem 6 Let I be a bounded closed interval of the X-axis which does not contain points Xj. n . Then 
there exists 7 = j(uij, I) > and a thin family of sets 1(e) such that the asymptotic expansion U5\) 
holds on all (finite and infinite) channels Cj. £ uniformly in X G I \ 1(e) and x in any bounded region 
of R d . The function ^ in Mty is a vector function on the limiting graph which satisfies equation 
conditions J17p at infinity, and the GC H8\) . 

Remarks. 1) It will be shown in the proof of Lemma [Til that for spider domains the estimate 
of the remainder is uniform for all x G R d . For general domains, we provide the estimate of the 
remainder only in bounded regions of R d in order not to complicate the exposition. 

2) The arguments, used to justify the asymptotic behavior of the scattering solutions and prove 
TheoremlHl can be applied to the study the asymptotic behavior of the outgoing solutions of the non- 
homogeneous problem ([5]) as e — > 0, A > Ao- The asymptotics will be expressed in terms of solutions 
of the corresponding non-homogeneous equation on the limiting graph. One can easily show that 
the GC can not be chosen independently of / even if we consider only functions / with compact 
support. However, if the support of / is separated from the junctions then the solution of the non- 
homogeneous equation on the limiting graph satisfies the same universal GC (|18[) that appear when 
scattering solutions are studied. The latter is related to the following fact: the outgoing solution in 
a narrow channel behaves as a combination of plane waves plus a term which decays exponentially 
outside of the support of / when e — > 0. 



6 



Note that the GC for the function on the limiting graph depend on A. In fact, there exists an 
effective matrix potential on T which is independent of A, and allows one to single out the scattering 
solutions 5 on T with the same scattering data as for the original problem in f2 E . These results will 
be published elsewhere. 

The convergence of the spectrum of the problem in £! e to the spectrum of a problem on the 
limiting graph has been extensively discussed in physical and mathematical literature (e.g., [I]-[7], 
O [El HSl [TBI HH] and references therein). What makes our paper different is the following: all the 
publications that we are aware of, are devoted to the convergence of the spectra (or resolvents) only 
in a small (in fact, shrinking with e — ► 0) neighborhood of Ao (bottom of the absolutely continuous 
spectrum), or below Ao. Usually, the Neumann BC on 9f2 e is assumed. We deal with asymptotic 
behavior of solutions of the scattering problem in Q. e when A is close to A > Ao , and the BC on dfl £ 
can be arbitrary. 

In particular, papers [5], [T2], [T3], [TS] contain the gluing conditions and the justification of the 
limiting procedure e — > near the bottom of the spectrum Ao under assumption that the Neumann 
BC is imposed at the boundary of Q e . Note that Ao = for the Neumann BC. Typically, the GC in 
this case are: the continuity of s (s) at each vertex v and 53y=i ( v ) = 0> i- e - the continuity of both 
the field and the flow. These GC are called Kirchhoff's GC. In the case when the shrinkage rate 
of the volume of the junction neighborhoods is lower than the one of the area of the cross-sections 
of the guides, more complex energy dependent or decoupling condition can arise (see [9], [13j . [7] 
for details). Let us stress again that this is the situation near the bottom Ao = of the absolutely 
continuous spectrum. As follows from Theorem [3) the GC and the small e asymptotics are different 
when A > Ao- 

Both assumptions (A — > Ao, and the fact that the BC is the Neumann condition) in the papers 
above are very essential. The Dirichlet Laplacian near the bottom of the absolutely continuous 
spectrum Ao > was studied in a recent paper [16 under the condition that the junctions are 
more narrow than the tubes. It is assumed there that the domain f2 e is bounded. Therefore, the 
spectrum of the operator (JTJ) is discrete. It is proved that the eigenvalues of the operator (fl]) in 
0(e 2 )-neighborhood of Ao behave asymptotically, when e — ► 0, as eigenvalues of the problem in the 
disconnected domain that one gets by omitting the junctions, separating the channels in fi e , and 
adding the Dirichlet conditions on the bottoms of the channels. This result indicates that the waves 
do not propagate through the narrow junctions when A is close to the bottom of the absolutely 
continuous spectrum. A similar result was obtained in [3] for the Schrodinger operator with a 
potential having a deep strict minimum on the graph, when the width of the walls shrinks to zero. 

We also studied the Dirichlet problem for general domains Q £ without special assumptions on 
the geometry of the junctions when, simultaneously, e — > 0, A — > Ao, and the diameters of the guides 
and junctions have the same order 0(e). Our conclusion is that, generically, waves do not propagate 
through the junctions when the frequency is close to the bottom of the absolutely continuous spec- 
trum. Let us stress that this is true both in the case when the diameters of the junctions are smaller 
than the diameters of the guides, and in the case when they are larger. Some special conditions 
must be satisfied for waves to propagate if A — * Ao. An infinite cylinder, which can be considered 
as two half-infinite tubes with the junction of the same shape, can be considered as an example of 
a domain where the propagation of waves at A = Ao is not suppressed. Less trivial examples will 
be given in our next paper. We do not deal with the problem near the bottom of the absolutely 
continuous spectrum in this publication. A detailed analysis of this problem will be published else- 
where. However, we show here that the GC on the limiting graph with A > Ao, generically, have a 
limit as e — > 0, A — > Ao, and the limiting conditions are the Dirichlet conditions. To be more exact, 
the following statement will be proved. 

Theorem 7 1) Assume that the resolvent h!2\l does not have a pole at k — VAo- Then the scattering 
matrix U0\) , defined for A > Ao, admits an analytic in z = yA — Ao extension to a neighborhood of 
the point z = and is equal to —I at z — 0, where I is the (mo x m$)-identity matrix and mo is the 
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number of infinite channels Cj. £ with Xjfl = Xq. 

2) Assume that the resolvent of the auxiliary problem in the spider type domain Cl VyS does not 
have a pole at k — yAo- Then the GC U8\) have limit as A — > Ao of the form 

eT' v j t ^\t)+2i^ v \t)=Q, t = 0, 

where T' v — -j^T v . The GC also have limit when e — * 0, A — > Ao independently. This limit is the 
Dirichlet condition ^ v \0) = 0. 

A simple version of the results presented in this paper (for models admitting the separation of 
variables) was published in our paper [14j . 

The next section contains the proofs of Theorem 2] and [31 The statements of these theorems 
mostly concern problems with a fixed value of e. Without loss of the generality, one can assume that 
e = 1 there. The last section is devoted to the proof of Theorem [S] on asymptotic behavior of the 
scattering solutions as e — > 0. Here the dependence of all objects on e is essential. At the end of the 
last section, one can find a proof and a short discussion of Theorem [7] 



2 Analytic properties of the resolvent R\. 

We denote by fie the following bounded part of fl £ : 

nw = n e \ u {c jJE n {* > a}). (19) 

j<m 

The next lemma will be needed later. 

Lemma 8 If the homogeneous problem @) with a real A > has a non-trivial solution u, then 
either \f\ is an eigenvalue of — e 2 A and u decays exponentially at infinity, or A 6 {Aj, n } and i U j 1 )) 
holds. 

Proof. From the Green formula for u and u in the domain fl £ a \ a > 0, it follows that 

Im f —TidS = 0. 

Jdn { E a) ov 

where v is the unit normal to d£l e °^ and dS is an element of the surface area. Using the boundary 
condition ([5]) we arrive at 



l / utudy = 0. (20) 



Im 

ldn { E a> \dfl e 

This, (J6]), and the orthogonality of the functions ipj >n imply, for a — > oo, 

V^- X j,n\a 3 ,n\ 2 + 0(e-^)=Q, 

j,n:\j jTl <\ 

which justifies the lemma after taking the limit as a — > oo. This completes the proof. 
Let C'a e be the channel Cj jS extended along the whole t axis, 

C' jj£ = {(t,ey):t€R, y G Uj C R n ~ 1 } 

We denote by the resolvent (fTTj) of the operator — e 2 A in the extended channel Cj . Let T\(C'^ ) 
be the set of functions from L 2 (C'j £ ) with the support in the region \t\ < a, and let H 2 (Cj e ) be the 
Sobolev space of functions in the domain C'j E n {b < \t\ < b + 1}. Consider the operator 

■■ L 2 com {C' 3 J - Ll c (C' j>£ ). (21) 
The following lemma can be easily proved using the method of separation of variables. 



Lemma 9 (1) The operator \21\) admits an analytic continuation from the upper half plane ImX > 
onto the real axis with the branch points at X — Xj. n , n = 0, 1, ... . 

(2) If Xj. mj < X < Aj >rnj +i and h G L^ oro (C'-) then R^h has the following behavior as t — > ±oo 

raj 

4 j) ^Ei e = If '^(^) + 0(e-*)''l), 7 >0, (22) 

n=l 

where 

4,n = = „■ y ! = / e^^^^MeMT, y)drdy. (23) 

*!y A - As," •'Wj.e oo 

($,) £e£ A £ (, where I is a bounded closed interval of the real axis such that Xj, mj < A < Xj^ mj +i 
for all A € I. Let h G Lf s (Cj E ) and b > 0. TTien t/iere exist positive constants c = c(l) and 7 = 7(2) 
which are independent of X G I, e and h, and such that the remainder term r in the right-hand side 
of H22]) has the estimate 

\\r\\ HHcU <ce-^\\h\\ LUc , E) . 

Proof of Theorem [4j The statements of the theorem mostly concern the problem with a fixed 
value of e. Without loss of generality, we can assume that e = 1, and we omit e in the notations of 
all objects (f2 e , C Jl£ , and so on). The dependence on e will be restored in some parts of the proof, 
when this dependence on e is essential. 

Step 1. Construction of the resolvent. Let us introduce the following partition of unity on fl 

m 

X> = 1. (24) 
3=0 

We fix arbitrary functions <pj G C°°(fi), 1 < j < m, such that <pj = 1 in the (infinite) channel 
Cj for t > 2, <p 3 ■ = in Cj for t < 1 and outside of Cj. The function <fio is defined as follows 
(f)Q = 1 — X)j<m "Ar We a ^ so nee d functions ipj that are equal to one on the supports of (pj, which 
will allow us to smoothly extend functions defined only on infinite channels or only in a bounded 
part of f2 onto the whole domain f2. We fix functions ipj G C°°(0), 1 < j < m, such that ipj = 1 in 
the infinite channel Cj for t > 1 (i. e. on the support of <pj), ipj = outside of Cj. Let ipo G C°°(fl) 
be a function such that ipo = 1 on the support of <f>o, and V'o = in all infinite channels Cj when 
t > 3. Note that 

^ij^j = <pj, < j < m. (25) 
We construct the parametrix (almost resolvent) for the problem ([5]) in the form 

m 

P A : L 2 (fi) -> L 2 (tt), P A / = V> O iM0o/)+X>i4%;/)- (26) 

where i?v is the resolvent (fTTj) of the operator in with a fixed A' = ia, a > 0, which will be chosen 
later, and are resolvents of the negative Dirichlet Laplacians in Cj. If / G L 2 (Q) then <pjf = 
outside Cj, and we consider <£j/ as an element of L 2 (Cj). Then the operator Rjp can be applied 
to (pjf and Ryi\(f>jf) G L 2 (Cj). Since ^ = at the bottom of Cj and outside of Cj, we consider 

ipjR^\(pj f) as an element of L 2 (il) that is equal to zero outside of Cj. In this way, the operator P\ 
is well defined for A ^ [0, 00). 

Let us look for a solution u G L 2 (Vl) of the problem (HJ with A ^ [0, 00) in the form of u = P\h 
with unknown h G L 2 (fl). Obviously, u satisfies the Dirichlet boundary condition since each term in 
(|26[) . applied to any h, satisfies the Dirichlet boundary condition. The substitution of P\h for u in 
equation ^ with A ^ [0, 00) (and e = 1) leads to 

(-A - A)P A /i = -(Aip )[Rx>(4> h)} - 2V^o • V[Ra'(M)] 



!) 



-^o(A + X')[R X , (cj> h)} - (A - \')MRx'^oh)} 

m 

- {(A^)[R^(^f)} + 2V^ • VR^fah) + ^-(A + XftBWfoh)]} = /> 

Using (|25[) . the last relation can be rewritten in the form 

h + F x h = f, (27) 

where 

F x h = -[(A + A-A')^][i^'(M)]-2V^-V[JiA'(M)] 

m 

- ^ {(A^)[^V^)] + 2V^ ■ V4 i} (^/i)} . (28) 

Let us show that the operator 

Fx '■ L 2 (fl) — > L 2 (il), A^[A ,oc), (29) 

is compact and depends analytically on A. Indeed, the resolvents i?A' and R^ map any function 
/ e L 2 into the solution of the problem (O in the domains f2, Cj, respectively. Thus, these operators 
are bounded as operators from L 2 into the Sobolev spaces H 2 . Since the formula (|2"5)) contains at 
most first derivatives of the resolvents, the operator F\, A ^ [A , oo), is bounded if it is considered 
as an operator from L 2 (f2) into the Sobolev space iJ 1 (il). Since V?/>o = V-0j = at points x E Cj 
with t > 3, from (|2"5|) it follows that, for any infinite channel Cj, 

F x h = 0, x € Cj R {* > 3}. (30) 

Hence, the Sobolev imbedding theorem implies that the operator ([29]) is compact. The analyticity 
of the operator ([29]) is obvious since the operators R^ depend analytically on A, and Rx> does not 
depend on A. 

Now we put A = A' = icr and show that ||-Fio-|| - > as a — > oo. In fact, since the norm of the 
resolvent does not exceed the inverse distance from the spectrum, we have that 

\\R$\\<l/v, (31) 

where the first norm is considered in the space L 2 (Q) and the second one is in the space L 2 (Cj). 
Multiplying the equation (O, considered in the domain Q or Cj, by u and integrating over the 

domain, we get the following relation for the functions u = Ryf and u = RV, f, respectively: 

||Vu|| 2 2 - icr ||ti||| 2 = J ufdx, 

which implies that 

||Vu||f. < | J ufdx\ < llulli.il/IU.. 

Thus, 

WRx'fWmm, II^/IIhmc,) < Ca-^WfWv. (32) 

Since the formula (f2"5| contains at most first derivatives of the resolvents, estimates (f3~Tj) . (|3"2"j) imply 
that ll-Fiffll — y as u — > oo. 

We fix A' = ia in (|26|) in such a way that ||-Fa'|| < 1- Then from the analytic Fredholm theorem 
it follows that the operator 

(E + Fx)- 1 :L 2 (n)^L 2 (n), A^[A ,oo), (33) 
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exists and depends meromorphically on A. From here, (|26[) and (|27[) the following representation for 
the resolvent follows 

R x = P X (E + F x )-\ A^[A ,oc). (34) 

Step 2. Analytic continuation of the resolvent. In order to extend the operator (|12p meromor- 
phically into the lower half plane Imfc < we need to repeat the arguments used to justify (j34f . 
Consider the space L 2 a {£l) of functions / € L 2 (Q) with supports in (see (HU), i.e. / = in the 
infinite channels Cj when t > a. Let / £ L„(f2). Without loss of generality, one can assume that 
a > 3. Then (|2"T|) and ([BIT]) imply that ft is also supported in fi( a \ i.e. F\ can be considered as an 
operator in L 2 (fl) : 

F A : -> L 3 a (Q), Ag[0,oo). 

Let x = Xo(i) be a function equal to one when t < a and zero when t > a. From Lemma [9] it follows 
that the operators 

X R[S : L 2 a (C 3 ) - L^C,-), Imfc > 0, 

admit an analytic continuation into the lower half plane with the branch points at k = ±y/Xj~n. 

Further, u = RjS f satisfies equation ((SJ) with A = k 2 for all complex k 6 C, and therefore the 
operators 

X R®, XVR$ ■■ L 2 a (C 3 ) - L*(C,), fc € C, 
are compact and analytic in the complex plane C. Since \ = 1 011 the supports of V"0j, < j < m, 
we can insert the factor x on the left of all the resolvents in Q28j) . From here it follows that the 
operator 

F k2 : L 2 a (Q) -» L2(fi), kC, 
is compact and analytic with branch points at k = ±^/Aj jn . Hence, the operator 

(E + F^)- 1 :L 2 a (n)^L 2 a (n), keC, (35) 

is meromorphic with the branch points at k = ± \f\j~n- Together with ()26|) , (I34|) and the analyticity 
of the operators i?^ : L 2 (Cj) — * Lf oc (Cj), k 6 C, this implies that the operator (fT2|) admits a 
meromorphic continuation to the lower half plane with the branch points at k — ±a/X^ and poles 
determined by the poles of the operator (|35|) . Obviously, the poles of the operator ([35]) may have a 
limiting point only at A = oo. 

Step 3. Spectral analysis. First of all note that the existence of the meromorphic extension of 
the operator (fl"2|) together with the Stone formula immediately imply that the operator H = — A 
does not have singular spectrum. The proof of this fact can be found in [17] (see Theorem XIII. 20). 

In order to prove the part of statement (1) of the theorem concerning the absolutely continuous 
spectrum of the operator H = —A, we split the domain Q, into pieces by introducing cuts along 
the bases t = of all infinite channels. We denote the new (not connected) domain by W , and 
denote the negative Dirichlet Laplacian in Q' by H' , i. e. H' is obtained from H by introducing 
additional Dirichlet boundary conditions on the cuts. Obviously, the operator H has the absolutely 
continuous spectrum described in statement (1) of the theorem. Thus, it remains to show that the 
wave operators for the couple H, H' exist and complete. The justification of the existence and 
completeness of the wave operators can be found in [T]. Another option is to derive the latter fact 
independently using Birman theorem stating that the validity of the inclusion 

(H - A)"" - {H' - X)- n e Ji (36) 

for some A and n > 1 implies the existence and completeness of the wave operators. Here J\ is the 
space of operators of the trace class. The inclusion (|36|) can be derived from (f34|) and a similar 
formula for the resolvent of the operator H' . This completes the proof of the statement about the 
absolutely continuous spectrum. 
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The discreteness of the set {A J ' e } of eigenvalues follows from the fact that the operator (TT2"]) is 
meromorphic in A and has poles at {X^ 5 }. The existence of the poles at {A 3,e } can be derived from 
the Stone formula. Another proof will be given below. 

Let us prove the part of statement (1) concerning the spider domains. If £l E is a spider domain, 
then there exists a point x(e) and an e- independent domain Q such that the transformation (see ([2])) 

L £ : x — > x(e) + ex, (37) 

maps into f2 e . In order to stress the fact that the operator H = — e 2 A in the domain f2 e depends 
on e, we shall denote it by H^ £ \ The operator —A in the domain Q shall be denoted by 
Obviously, 

= L £ H^LJ\ (38) 

and this implies the independence of the eigenvalues of the operator of e. This completes the 
proof of statement (1). 

Step 4, real poles of the resolvent. The first part of statement (2) about the existence of the 
analytic extension of the resolvent was justified in step 2 of the proof. Now we are going to prove 
the second part of that statement concerning the set of real poles of the operator (fT2"|) . We denote 
this set of poles by K. Let us assume that either u is an eigenfunction of the operator H = — A 
with an eigenvalue A = A' > or u is a non-trivial solution of the homogeneous problem ([5|) , (fT3|) 
with A = A' > (recall that we assume that e = 1). We are going to show that k = zL^/xJ £ K. 
Consider the restrictions Uj of u to the cylinders Cj, 1 < j < m. Let Vj £ L 2 (Cj) be the solution of 
the problem 

(—A — X)Vj =0, x £ Cj] Vj = on d Cj, Vj = Uj when t = 0, 

where A ^ [0, oo), and d'Cj is the lateral boundary of Cj. The solution Vj £ L?{Cj) of this problem 
is unique and can be found by separation of variables. The function Uj satisfies the same equation 
with the fixed A = A' and the same boundary conditions. It is also defined uniquely by its values at 
t = and can be found by separation of variables. This implies that Vj converges to u as A — > A' + iO. 
Since u is a solution of a homogeneous elliptic problem, u £ C°° . Thus, Uj is infinitely smooth when 
t = 0, and the convergence Vj — > Uj takes place, for example, in the Sobolev space H 2 on the part 
of the cylinder Cj where < t < 2. Let 

m 

v = Y^ <§>i v i + ^ou £ L 2 (Cj), A £ [0, oo), 
j=i 

where {4>j} is the partition of unity which was introduced above. The function u can not be equal 
to zero identically on f2\ U Cj due to the uniqueness of the solution of the Cauchy problem for the 
operator —A — A'. Thus 

\\vWmn\uCj) = |MU 2 (n\uQ) = co > 0. (39) 

On the other hand, 

m 

(-A -X)v = - ^2[(A(j)j)vj + 2V(j)j ■ Vvj] - (A - \')<j> u - (A0 O )« - 2V0 O • Vit. 

Thus, (—A — X)v £ L 2 (fl). From the convergence Vj — > u and (|24|) it follows that (—A — X)v tends 
to zero in L 2 a {Q) as A — > A' +i0. Together with ([3^)1 this provides the existence of the pole of the 
operator (|12[) at k = y/X' '■ The pole at k = —VX* exists due to the relation R\ = Rj. Hence, 



12 



Now let us assume that at least one of the points ±y/X/ belongs to K. The relation R\ = Rj 
implies that the second point also belongs to K, i.e. ^f~X! £ K, and there exist a > and / S 
such that 

W;=flA/= (A^ ^A^r-' ; IHI^ ( o(^))< C , A^A' + zO, (40) 

where n > 1 and u does not vanish identically. In fact, n can not exceed one, but it is not important 
for us now. Obviously, 

(-A-A')« = 0, x£fl; « = onffl. (41) 

From here and Lemma [8] it follows that in order to complete the proof of the second statement of 
the theorem it is sufficient to show that the asymptotic expansion ^ holds for the function u. 

Note that (|41 [) implies that u S C°° . Since / = in all infinite channels Cj when t > a, from 
relation ([30]) it follows that 

(-A - X)v = (A - X')u, xeCjD{t>a}; v = on dn. 



From here, the estimate in (|40j) , and standard local a priory estimates for solutions of elliptic 
problems it follows that for any vector a 



||^|<c(a), xeC j n{a+^>t>a+^ 

and therefore 



d a [(\-x') n w] ^ d^u 

dx a dx a [ ' 

uniformly on Cj fl {t = a + 1} as A — ► A' + iO. We restrict the functions (A — X') n w and u to 
Cj n {t = a + 1} and expand the restrictions with respect to the basis {<ySj>} of the operator —A 
in the cross section of the channel Cj. Let 7j, n (A) and 7° n be the coeficients of these expansions. 
Then H2]) implies that for any (3, 



|7^(A)-7°„I <cpn-P, X^X' + iO. (43) 
The function w := (A — X') n w satisfies the following relations in Cj D {t > a + 1} : 

(-A-A)w = 0, w = for x £ dCj n {<> a + 1}, w| t=0+ i = lj,nW<Pj,n(y), 

n 

where A <^ [0, oo). One can find the solution w £ L 2 of this problem by the method of separation 
of variables and then pass to the limit as A — * A' + iO using (|43|) . This leads to the asymptotic 
expansion ([6]) for u and completes the proof of the second statement of Theorem 01 

Step 5, the proof of the last two statements of the theorem. If k = V^Y, A' > 0, is not a pole or a 
branch point of then 

w:=R x f = u(x) + (X-X')v(x > X); \ \v\ | £a(tl c.)) < c(o), A -> A' + i0, 

where a > is arbitrary and 

(-A - X')u = f, igO; u = on 90. 

In order to prove the third statement of the theorem, we need only to show that the asymptotic 
expansion ([6]) holds for u. It can be done exactly in the same way as it was done for function u in 
(HO]) by representing u in Cj n {t > a + 1} as the limit of functions ic as A — > A' + iO. 
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In order to prove the last statement of the theorem one can look for the solution '5 = ^ s ^ of 
the scattering problem in the form 



* = cj> s e- i ^ x - x >' kt <p Sik (y)+u, 

where 4> s is the function from the partition of unity This reduces problem |(7J), © to the 

uniquely solvable problem (O, ([6]) for u. 
This completes the proof of Theorem [U 

Proposition 10 Let <J s ^ and ^ s ',k' be two scattering solutions, and let a*'^ be the transmission 
coefficients for the scattering solution ^ St k- Then 

1) The following energy conservation law is valid: 

m m i 

2) If these solutions correspond to different incident waves ((s,k) ^ (s',k')), then 

Proof. Since the statement concerns the problem with a fixed value of e, one can put e = 1 and 
omit e in the notations O g , fig . Green's formula for ^s^. and , 3>s',fc' in the domain fiW implies, 
similarly to (p?0|) , that 

/ [(* s .fc)t* S ',fc' - * s ,fc(* s ',fc')t] d i/ = o. 

From here, ((5J), and the orthogonality of the functions ^ „ it follows that 



En \ s,k~s',k' f\ T — s,k — 2i»/A — A« fe a 

V A - V°i, n a j,i _ V A ~ X ^ a j,n e V 

J," 

+ \/A- \ a >,k>a.*'*' e 2 V A -V. fc '« _ - A s , fe £ + 0( e -^ a ) =0, a ^ oo, 

where 5 = 1 if ^ ^ ^ = f ^/ ^ ) an( i <5 = otherwise. We take the average with respect to 
a e (A, 2A) and pass to the limit as A — > oo. Then we get 



which justifies both statements of the proposition. This completes the proof. 

Proof of Theorem^ PropositionfTUlis equivalent to the relation A* A = I for A = D 1 ' 2 TD- 1 / 2 , 
which provides the unitarity of the matrix A. If one applies Green's formula to the scattering solu- 
tions $ Si j; and ^ s ',k' i then the arguments used in the proof of Proposition [TU] lead to the symmetry 
of D 1 I 2 TD- 1 / 2 . This completes the proof of Theorem O 

3 Asymptotic behavior of scattering solutions as e — > 0. 

We start with a study of scattering solutions in spider domains f2 e 
Lemma 11 Theorem [3| is valid for spider domains. 
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Proof. The transformation L~ l , see (I57|) . maps the spider domain fl £ into the e— independent 
domain with the channels Cj, 1 < j < m, The coordinates (£, y) in Cj are related to coordinates 
(t, y) in C^e via the formulas 

t = t/s, y = y/e. (44) 
The scattering solution ^ = $ s of the problem in 17 has the form similar to ©: 

n=0 

x e Cj, t —y oo. 

Since ^ s ,k is a smooth function, the remainder term r in the formula above can be estimated for all 
values of t : 

\r\ < Ce~^, xeCj. 

Since the scattering solutions in the domains £l e and are related via the formula ^ s- k{x) = 
1 $> St k(Lj 1 x), it follows that 

*S = ^ J e _l ^ ii V s ,fc(y/e) + ^^,r l e l ^^ i V J ,r l (y/e) + r (e) , 

n = 

|r( £ )| < C*e- 7 */ e , xeCj. (45) 

Thus, the asymptotic expansion p3|) , (JTTJ) is valid, and it only remains to show that the GC (flB]) 
holds for vectors c = <; s .k determined by (|45[) (the definition of these vectors is given in the paragraph 
above formula (|18|) V We form the matrix £ = £(t) with columns <; s _j. taking them in the same order 
as the order chosen for elements in each of these vectors (first we put columns with s = 1 and 
A; = 1,2,..., mi, then columns with s = 2, and so on). From (|4"5"1) it follows that 

E(0)=J + T, S'(0) = iD(_J + T), 

where T is the scattering matrix, / is the identity matrix of the same size, and D is the diagonal 
matrix of the same size with elements \j, n on the diagonal. Hence, e(I + T)D _1 S'(0) + — 

T)E(O) = and GC (jT5]) holds for the columns of the matrix S. This completes the proof of the 
lemma. 

The following two lemmas about spider domains will be needed in order to prove Theorem [6] for 
general domains. Let R\ be the resolvent of the operator H = — e 2 A in a spider domain f2 E , and let 
R\ be the resolvent of the similar operator H = — A in the domain £1 which is the image of fl £ under 
the map Lj 1 , see (|37[) . Note that the operator R\ and its domain, L 2 (£l £ ), depend on e, while the 
operator R\ is e- independent. Formula (|38p implies 

Lemma 12 The following relation holds 

R\ = L e R\L s , 

Let us fix m constants tj > 0, 1 < j < m. Let Vt e be a spider domain with the channels Cj jE , 
1 < j < m. Consider slices D^ e of Cj <e defined by the inequalities |t — tj\ < 3e. Let il' £ be a bounded 
domain which is obtained from f2 £ by cutting off the infinite parts of channels Cj. £ on which t > 
Let a function h G L 2 (fl e ) be supported in one of the domains -Dj, £ , for example, with j = s. 

Below, when the resolvent R\ of the operator H = — e 2 A in Sl £ is considered with A belonging 
to the continuous spectrum of the operator, R\ is understood in the sense of the analytic extension 
described in Theorem [4] We denote the Sobolev spaces of functions which are square integrable 
together with their derivatives of up to the second order by H 2 (fl' £ ) and H 2 (Dj iE ) . 
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Lemma 13 Let fl E be a spider domain. Let I be a bounded closed interval of the X-axis that does 
not contain points Aj jTl , and let a function h G L 2 (fl £ ) be supported in the domain D s e . Then there 
exists 7 = 7(wj, I) > such that 

(1) Rxh = J2 c *-k*s,k+r in tf e , |r (x)|<-— — —\\h\\ LH n e) , (46) 

where "J^fc are scattering solutions, the coefficients c s ,k — c~ k (h) are given by i23\) , and \ J are 
eigenvalues of the operator H in Q, e (see statement (1) of Theorem^: 

m s rrij 

(2) /,v, /, . X'-^XC'X ~ " ' ■ "> l} - ■ ( 4? ) 

fe=0 ra=0 

where \ \r 3 \ \ H 2 (D . e) < n | A _ ^jj INU 2 (Q e )- 

Here S s j is the Kronecker symbol (6 s .j = 1 if s = j, S s j — if s ^ j), is the resolvent of —A in 
the extended channel C' s e (channel C s>e extended to — oo along the t axis), c s .k — c~ k (h), and t s -' n 
are the transmission coefficients (see the remark following definition^. 

Proof. Let a function a € C°°(fl s ) have the form: a = 1 in C SyE when t > |ij + e, a — in 
f2 e \C S)£ , and a = in C S)£ when t < ^tj. Consider the function 



u = aR { ° ] h + c s . fc [* s ,fc - ae _i — H*-V.,fc(v/e)]. (48) 

fe=0 

Obviously, u — on dfl £ , since each term in the right hand side above satisfies the Dirichlet boundary 
condition. Furthermore, 

- e 2 Au - Xu = h-e 2 [Va- VR^h + (Aa)R { ° } h - Va • \7g - (Aa)g\, (49) 

where 

fe=0 

The right hand side in (|49[) has the form h + hi, where hi is supported in the slice |i s < t < |t s + e 
of C Sj£ . From Lemma [9] it follows that 

\\hi\\ L z ( n e )<Ce-^\\h\\ L z ( n e ). 

It is also clear that the behavior of the function u at infinity is described by j6|. Hence, u = 
R\(h + hi) due to the statement (3) of Theorem |U From here and (|48|) it follows that 

R x h = aR[ s) h + Cs.k(*s,k ~ ae-'^^^s^y/e)) - R x h x . (50) 

This implies equality (|46|) with ro = —R\hi. Let fi" be obtained from J7 e by cutting off the parts 
of channels Cj >e where t > ^tj. Since operator (fT2")l is meromorphic (due to Theorem 2]) and has 
poles of at most first order due to the Stone formula, | |^"o | |i 2 (n") can be estimated by the right 
hand side of inequality (|46|) . Since (— e 2 A — A)ro = in £1" and ro = on the lateral side of <9£1", 
standard a priory estimates for elliptic equations lead to the estimates on ro in Sobolev norms in Q' E . 
These estimates, together with Sobolev imbedding theorems, justify the estimate (|46l) . Similarly, 
(|47| follows from ([50]) and Lemma fT2l This completes the proof of the lemma. 
We need two more auxiliary statements in order to prove Theorem [5] 
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Lemma 14 Let a real-valued function f belong to C n+ (R ) and 

n+l 

| C „ +1 =^sup|/W|=^ + <oo, (51) 



k=0 



\f {k) (x)\ > A- > 0, xeR 1 . (52) 

Then for any a < A~ /2, i/ie set T a = {x : |/(x)| < cr} has the following structure. There exists 
a constant c which depends only on A ± and n and such that, for any bounded interval A C R 1 , 

a) the number of connected components ofT a in A is finite and does not exceed c(|A| + 1), 

b) the measure of each connected component ofT a in A does not exceed ca x l n . 

Remark. The last estimate can not be improved. In fact, if f(x) = sin" x then T a [— — ? ] 
- 2a 1 /™. 

Proof. We shall denote by Cj different constants which depend on A ± and n but not on /. If 
x € T a then (|52|) implies that 



J2\fW{x)\>A-/2, 

k=l 

and therefore, |/^ fc '(a;)| > A~/2n for the chosen x and some k = k(x), 1 < k < n. Since |/*- fe+1 - ) | < 
A + , x £ R 1 , there exists an interval A x such that x £ A x , |/' fe )(ir)| > A~/4n on A x , and | | — 
Co = 4 ^ A + ■ The set of intervals A^ covers T a n A. Hence, one can select a finite number of intervals 
A x covering r„n A. Then one can omit some of them in such a way that the remaining intervals 
still cover r CT H A with multiplicity at most two. This leaves us with at most + 1) < Ci(| A| + 1) 
intervals A^ covering r CT n A. Thus, it is enough to prove the lemma for an individual interval A' 
(one of the intervals A x ) such that |A'| = Co and \f^ k '(x)\ > C2 on A' for some fixed value of k, 
1 < k < n. 

Equations f(x) = ±cr have at most k solutions on A'. In fact, if there exist k + 1 points where 
f{x) — a then there are k intermediate points where f'(x) = 0. Thus, there are k — 1 points where 
f"(x) = 0, and so on. Finally, there has to be a point where f( k >(x) — 0. This contradicts the 
assumption that > C2 on A'. Hence, the set T a n A' consists of at most k + 1 intervals. It 

remains only to show that the length of these intervals does not exceed ca x l k . In order to estimate 
this length, we assume that there is an interval \x\,x\ + h] where \f{x)\ < cr, \f (k) {x)\ > c 2 . Put 
h' = h/k and consider the fc-th difference 

A fc = f( Xl ) -(f) f(*x + h') + ( * ) f( x i + 2h ') - - + (- 1 ) fc /( a; i + kh ')- ( 53 ) 

There exists a point £_ k e [x 1 ,x 1 + h] such that A fc = (h') k f<- k > (£ k ). Thus, |A fc | > ^ = c 3 h k . On 
the other hand, from (155)) and the estimate \f(x)\ < a it follows that \A k \ < a2 k . Hence, c^h k < a2 k , 
i.e. h < ca 1 ^ . This completes the proof of the lemma. 

Lemma 15 Let a set of functions f e — / e (A), e — > 0, on a closed interval I C R , have the form 

M 

fe = ^2Cj(X)e 1 ^- , (54) 
3=1 

where functions Cj(X) and real valued, functions gj(X) are analytic, there are no two functions gj{\) 
whose difference is a constant, and 

M 

£|C,(A)|>1. (55) 

3=1 
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Then, for any rj > 0, the set 

r v (e) = {\:\f e (\)\<e-^} (56) 
is thin (see the definition in the introduction). 

Proof. Consider the set Fo where g'i(X) = <?j(A) for some i ^ j. Due to the analyticity of the 
functions <?j(A), this set consists of a finite number points. Let us denote the number of points in 
T by a. Let V s be the <5/2-neighborhood of L°. Then 

\g' i (X)-g' j (X)\>a(5)>0, i^j, \el\T 5 . (57) 

Consider the functions 

3=1 

For any k, 

d k M < ) 

j-pUtJ) = ^^Mf^Me''^ + 0(e), (58) 

3=1 

We move the remainders to the left hand side and consider |58|) with 1 < A; < M as equations for 

unknowns Cj(e\i)e % s . The matrix of this system of equations with the elements Oj,fe = [gj[£fx)\ 
is a Vandermond matrix, and its determinant is bounded from below due to (|57| . This and ([55]) 
imply that 

j2\^m\>A-(s)>o 

3=1 P 

if e is small enough. It also follows from (|58|) that 

3=1 P 

Hence, Lemma [Til is applicable to at least one of the functions Re/ e (/x) or Im/ e (/i) on each connected 
interval of the set I \T S stretched by a factor of e -1 . Since we have at most ci + 1 those intervals, 
this implies that the set {A : |/e(A)| < a} can be covered by V s and C2(<5)e _1 intervals of length 
C2{5)a 1 ^ M . We take a — e~ v ^ e , and this completes the proof of the lemma. 

Proof of Theorem [6] The proof is based on a representation of the resolvent R\ of the problem 
in fl e through the resolvents R v .\ of the operators H = — e 2 A in the spider domains £l Vt£ , formed 
by an individual junction, which corresponds to a vertex v, and all the channels with an end at 
this junction, where the channels are extended to infinity if they have finite length. Let us consider 
the slices Dj E of the finite channels Cj. £ , j > m, defined by the conditions tj < t < tj + 3e where 
tj = Alj/b. We construct the following partition of the unity on £1 £ : 

X] ^ = 1 ' 

where V is the set of all the vertices v of the limiting graph r, (j> v S C°°(f2 e ), and is defined as 
follows. The function <j) v is equal to one on the junction J„, which corresponds to the vertex v, on 
the infinite channels adjacent to J v and on the parts of the finite channels adjacent to J v where 
t < tj + e. The function <p v is equal to zero on the parts of finite channels adjacent to J v where 
t > tj + 2e, and also on all the other junctions and channels which are not adjacent to J v . Let 
tl) v G C°°(0 £ ), tp v = 1 on the support of tp v = on the parts of finite channels adjacent to J v 
where t > tj + 3e, and also on all other junctions and channels which are not adjacent to J v . 
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We fix a vertex v = v. Let Jg be the corresponding junction of f2 e . We choose the parametrization 
on r in such a way that the value t = on all the edges adjacent to v corresponds to v. The origin 
(f = 0) on all the other edges can be chosen at any end of the edge. We are going to justify the 
asymptotic expansion (|15[) in the domain C(v) consisting of the infinite channels adjacent to Jy and 
the parts t < 3/j/5 of the finite channels Cj. £ adjacent to Jg. Moreover, it will be shown that the 
function <; in the asymptotic expansion satisfies equation (|16p . conditions (|17p at infinity, and the 
GC (|18p . Since v is arbitrary and the union of all domains C(v), v G V, covers all the channels, the 
validity of (fT5| in C(v) justifies the statements of Theorem [6] 

Let us show that the asymptotic expansion (fT5|) in C(v) for any scattering solution follows 
from a similar expansion for functions of the form u = R\f, where / G L 2 (fl e ) is supported in 
UDj^. In fact, let u — ip Vl ^s\ v x > wnere the vertex v\ — v\{s) corresponds to the first junction 

J Vl encountered by the incident wave, v is the solution of the scattering problem in the spider 
domain Cl vltS , and the function u is considered as a function in f2 e which is equal to zero outside of 
the support of ib vi . Then 

(s 2 A - X)u = /, / := -s 2 [^ Vl • V¥ s % tvi + (A^J*S, V J. 

Obviously, / G L 2 (f2 e ) and / is supported in UDj tS . From statement (3) of Theorem [4] it follows that 
there exists the unique outgoing solution v = R\f of the equation (— e 2 A — X)v = /, A G I \ {A- 3 }. 
Then 

since this function satisfies and © . From here and Lemma [TT] it follows that the asymptotic 
expansion ([15)) and the properties of <j mentioned in Theorem [S] hold for in C(v) if the corre- 
sponding properties are valid for R\f in C(v). Hence, the proof of the theorem will be complete as 
soon as we show that, for any / G L 2 (fl e ) with the support in UDj j£ , the function u = R\f has 
expansion <|T~5|> in C(v) with j3j >n = and ^ satisfying the GC (fl8|) . 
Consider the operator P\ defined by the formula 

P x h = ^ lM2«,A(M)i A e I, (59) 

where /i G L 2 (fi e ) is supported in Ufl 3i£ , Z is defined in the statement of Theorem [6l and the 
resolvents R v ,\ for real A G I are understood in the sense of Theorem |U We look for u — R\f in 
the form of P\h with an unknown h G L 2 (Q, e ). This leads to the equation (compare with (|27[) . (|28p) 

h + F x h = f, F x h=-s 2 J2[ 2 ^v^RvA<Pvh) + (Aip v )R ViX (cj> v h)]. (60) 

From here, similarly to (J34]) , it follows that 

Rxf = Px(I + F x y 1 f (61) 

for ImA > 0. Similarly to (f35|) . one can show that the operator 

(I + Fx)- 1 : L 2 a (n e ) ^ L 2 a (Q £ ) (62) 

admits a meromorphic extension into the lower half plane ImA < with the branch points at 
A = Xj .„. The only difference is that now we use operators R Vi \ instead of R x , and R v .\ depend 
meromorphically on A, while IvP are analytic in A. So, one needs to refer to a version of an analytic 
Fredholm theorem where the operator may have poles (with residues of finite ranks). This version 
of the theorem can be found in [5] , and applications of this theorem to operators similar to (f6"2"|) can 
be found in [T5], [2D]. Hence, formula (|6"Tj) is established for all complex A. 
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All the operators in (|6ip depend on e. The function (/ + Fx)^ 1 f is meromorphic in A, and its 
poles depend on e. In order to find a set on the interval I where the operator (/ + F\) exists 
and is bounded uniformly in A we shall use the following reduction of equation (|60p to a system of 
equations where the domains of the operators do not depend on e. 

Recall that / is supported in UDj t£ . Formula (|60| for F\ implies that the function F\h is also 
supported in UDj tS . Thus, the support of any solution h of (|60|) belongs to VJDj tS . We shall identify 
functions / and h with vector functions whose components are the restrictions of / and h, respec- 
tively, to individual domains -Dj, e , m + 1 < j < N. Furthermore, we map Dj e onto e-independent 
domain Dj by the transformation L' e defined by formulas 

t-tj = et, y = ey. (63) 

This transformation differs from (|4"4")l by a shift in t (compare with ([3T)) V The vector functions 
of variables (t, y) with components from L 2 (Dj) defined by /, h will be denoted by /' and h', 
respectively. Then equation (|60[) can be written in the form 



b! + F' x h' = /', 

where F' x is the \(N — m) x (N — m)]-matrix operator which corresponds to the operator F\. Here 
N — m is the number of finite channels in Vl e . Recall that the entries of the vectors /' and h' are 
functions with the domains Dj, which do not depend on s (and A). Our next goal is to describe how 
the entries of the matrix F' x depend on s and A. It will be done using (|60p . where each resolvent 
R v ^x can be specified using (|T7|) . 

The first term in the right hand side of (|47p. 



i?i s) :L 2 (C; e )-L 2 (C; e ), 
depends on e. The transformation (|63[) maps C' s e onto the e- independent cylinder C' s . The operator 

R[ s) :=L' e R^(L' e )-':L 2 (C' s )^L 2 (C' s ) 

does not depend on e (Lemma I12|) . and it depends meromorphically on A. Thus, the contributions 
from the first term in the right hand side of (|47|) to the entries of the matrix F' x are operators which 
are e-independent and meromorphic in A. The rest of the terms in the right hand side of (|47| (other 
than the remainder) are operators of finite ranks. Due to Lemma IT51 (see also the formula (I2U1) for 
Cs,k = c 7fc)> the contributions of these terms to the entries of F^ are e-independent operators which 
are analytic in A and are of the rank one, multiplied by functions q v -j, n ,s.k of the form 



. o 3v /A-A : ,- i „+/i sv /A-A s k 

3^;i,n,s,fc(A,e) = e % •= . (64) 

Here ay = tj or ay = lj — tj (independently, (3 S — t s or {3 S — l s — t s ). Formula (|47|) leads to a.j = tj, 
[3 S = t s if 1) the channels Cj i£ and C Si£ are adjacent to a common junction, which corresponds to the 
vertex v, and 2) the parameter t on both channels Cj jS and C Sj£ is introduced in such a way that t = 
at the vertex v. Other options in the choice of a.j and (3 S correspond to opposite parametrization of 
the channels Cj jS , C Sj6 , or both. If Cj^ e and C Si£ do not have a common junction which corresponds 
to the vertex v then q v; j t n tSt k = 0- Thus, the matrix operator F' x can be represented in the form 



^2 qv-j,n,s,k{\ £ ) F \ 



v.n.k 



+ R, (65) 

j,s>m 



where F®, Fi' n,s ' are e-independent operators, F® is meromorphic in A, operators F(' n ' s ' are 
analytic in A and have rank one, the summation extends over all the verticis v and over n € [0, rrij], 
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k G [0, m s ]. The operator R = R(e, A) corresponds to the remainder term in (jTT)) and has the 
following estimate 

||p||<_^L_ 

11 11 - n V6 ,|A-A»T 

Since the analytic Fredholm theorem |2j is applicable to the operator J + i 7 ^, from (f65|) it follows 
that it is also applicable to the operator I + F?. Let be the <5/2-neighborhood of the set consisting 
of both the poles A J of the operator (I + F®) -1 located inside / and the points A J £ !. Then 

Wil + F^W^CiS), \\R'\\<C(6)e-^, \el\l 5 , (66) 

where 

R' = R(I + F^y 1 . 

Formula (f65|) implies, for A G I \ I s , 

(I + F^)- 1 = (7 + F")- 1 [I + qG + R']- 1 , (67) 

where qG is the matrix operator with matrix elements ^2 v . n j, 1v,j,n,s,k{X, £)G^ n ' s ' , N — m < j,s < 
iV. Here 

The operators G^ n ' s ' k are meromorphic in A and have rank one. 
The equation 

(7 + gG)z - 5 (68) 

for x can be reduced to an equation in the finite dimensional space 5* spanned by the ranges of the 
operators G^ n ' s ' k . We fix a basis in S, reduce equation (|68|) to the algebraic system Ax = g for 
coordinates of the projection of x on S, and solve the system using the Kramer rule. Since functions 
Qv-j,n,s,k are bounded when A G I, the procedure described above allows us to estimate the norm of 
the operator (I + qG)^ 1 through | det -1 A\. 

Hence, there exist a polynomial P = P(q v -j, n ,s,k) of variables q v -j, n ,s,k which has the following 
properties. Its coefficients are meromorphic in A (with poles belonging to the set {A J } U {A J }), the 
polynomial is linear with respect to each variable q v -j.n,s,k , and is such that 

11(7 + qG)- 1 \\ < C\f e (\)\-\ Xel\l s 7 / e (A) := l + P(q v . ij>n , s , k (X,e)). 

The function / e (A) here has the form ([54]) with one of gj identically equal to zero, and the corre- 
sponding coefficient Cj equal to one. The latter implies (|55|) . Thus, Lemma 1151 can be applied to 
the function / e (A) above on each connected interval of the set I \ I s . There are only finitely many 
such intervals. Thus, for any r\ > 0, there exists a thin set T v (e) such that 

Wil + qGy'W^Ce^", Xel\T v (e). 

We choose r\ < 7, where 7 is defined in (|66[) . Then 

Wil + qG + R'y 1 ]] <Ce' /2e , \el\ T v (e), 

when e is small enough,. A similar estimates holds for operator (f67j) : 

||(/ + ^)- 1 ||<C7 e ^, AG/\r„( £ ). 

Hence, the same estimate is valid for the operator (I + Fx) , and from (f59|) . ([61]) it follows that 
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where h £ L 2 (fl e ) is supported in \JDj tS , j > m, and 

\\h\\ L ^ e) <Ce^ 2s \\f\\ LHne) , \el\Tr,(e). (70) 

Relations (|70p . (J3H1), an d Lemma [TT1 together provide the asymptotic expansion (TT5)) for i? A / 

needed to complete the proof of Theorem [Sj 

The last result, which we are going to discuss now, concerns the limiting behavior of the GC as 
A approaches Ao , the bottom of the absolutely continuous spectrum. We assume that the resolvent 
(fT2| does not have a pole at k = y/X~o. Obviously this assumption holds for generic domains f2 e . 
Theorem |4] implies that this assumption is equivalent to the absence of bounded solutions of the 
homogeneous problem (|5|) with A = Ao- Recall that the scattering matrix (flTJ| depends on A > Ao 
and the GC (|18[) depend on both A > Ao and e > 0. 

Proof of Theorem [7j Consider an infinite channel C s , for which A Si o = Ao- Let ^ s0 ^ e ^ ne 
scattering solution which corresponds to the incident wave 

ipinc = e 1 e (p a ,o(y/e). 

Let 4> s £ C°°(fl £ ), 4>s = 1 in the channel C s for t > 2, <p s = in C s for t < 1 and outside of C s . We 
represent in the form 

*i!o = 4>s^inc +U, A > A . 

Then u is the outgoing solution of the problem 

(-e 2 A -X)u = f, x € n e ; u = Q on dCl s , 

where / = — e 2 (A(j> s )ipi nc — 2e 2 Wcf>s Vipinc has a compact support. Hence, u = R\f. From here, the 
second statement of Theorem 01 and the absence of a pole at k = ^/Xq it follows that u, if considered 
as an element of L 2 oc (Vl e ), is analytic in z — y'A — Ao in a neighborhood of the point z = 0. Then 
from standard local a priory estimates for solutions of elliptic problems it follows that u is analytic, 
if considered as an element of any Sobolev space of functions on any bounded part of fl E . Hence, 
the restrictions Uj of u to cross-sections t = 2 of infinite channels Cj are analytic in z. Thus, for any 
infinite channel Cj , u is an outgoing solution of the problem 

(-e 2 A - X)u = 0, x e Q. e H {t > 2}; u = on9fi E fl{i>2}; u\ t =2=Uj. (71) 

Since Uj is analytic in z = yA — Ao in a neighborhood of the point z = 0, the coefficients aj.„ in the 
asymptotic expansion ([6]) for the solution u of (|71|) are analytic in z. This proves the analyticity of 
the scattering matrix. 

From analyticity of u in z and (fTTj) it also follows that the scattering solution when z = 0, 

is a solution of the homogeneous problem ((5]) with A = Ao, and satisfies (fT3|) . Thus ^ = when 
z = due to Theorem [U This implies that T — —I and completes the proof of the first statement. 
The second statement of the theorem is an obvious consequence of the analyticity of T v and (| 18|) . 
This completes the proof of the theorem. 

Remarks concerning Theorem [3 1) Consider a bounded domain il E with one junction and 
several channels of finite length. Let Q' s be a spider type domain which one gets by extending the 
channels of fl e to infinity. The spectrum of the problem |JSJ| in Cl e is discrete, and there exists a 
sequence of eigenvalues which approach Ao as e — > 0. Each of these eigenvalues has the form 

X n (e) = X + O(e 2 ). (72) 

Theorem [ll concerning the problem in W e , can be used to specify the asymptotic behavior (j72|) of 
the eigenvalues A„(e). The last statement of the theorem and (|72[) indicate that, for generic domains 
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f2 e , the asymptotic behavior of X n (e) as e — > (when n is fixed or n — > oo not very fast) is the same 
as for eigenvalues of the corresponding Dirichlet problem on the limiting graph with the Dirichlet 
GC at the vertex. This result will be discussed in more detail elsewhere. 

2) Our paper P3] contains a mistake in the statement of Theorem 5.1 (which is a simplified 
version of Theorem [7] above) about the form of the GC at the bottom of the absolutely continuous 
spectrum: k — > has to be replaced there by k — > 0, e — > 0. The arguments in the last 5 lines of the 
proof are wrong, but can be easily corrected with the additional assumption that e — > 0. (Also, the 
index of summation in the formulas (5.2), (5.4), (5.6) of that paper must be n, not j). 
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